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We compute the r-matrix for the elliptic Euler-Calogero-Moser model. In the trigonometric limit we 
show that the model possesses an exact Yangian symmetry. 
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1 Introduction 

The Euler-Calogero-Moser model was defined in p, 0] . In we considered the rational case and we derived 
the r-matrix. In this paper we are interested in its trigonometric and elliptic generalizations. In the elliptic 
case we compute the r-matrix and show that the usual elliptic Calogero-Moser model and its r-matrix are 
obtained by Hamiltonian reduction. In the trigonometric case we show that the current algebra symmetry 
exhibited by Gibbons and Hermsen p] in the rational case, is deformed into a Yangian symmetry algebra. 

We consider a system of N particles on a line with pairwise interactions. The degrees of freedom consist 
of the positions and momenta {jpi,qi)i=i---N and of antisymmetric additional variables (/^ = —fji)i,j=i---N, 
with the Poisson brackets 

{Pt,qj} = S,j (1) 

{fij, fki} — -^ [Sii f'jk + Ski fij + 5jk fii + Sij fki)- (2) 

The Poisson brackets of the fij just reproduce the 0{N) Lie algebra. The Hamiltonian will be taken of the 

form 

1 ^ 1 ^ 

H=^'^pi~^^Iio Io^V{qiJ), Qtj^qi-qj (3) 

i=l i,] = l 

with an even potential V{—x) — V{x). 

The equations of motion are easily derived: 

qt = Vi 

N 



Pi 



3 = 1 

N 



k = \ 
k^i.j 



Such a system admits a Lax representation only for specific potentials. Indeed writing the following ansatz 
for the Lax pair 



N N 



^(A) = ^Pi ^ii + X! ^('J'U' A) fij Cy (4) 

i=l ».3 = 1 

N 

M(A) = Em(9y,A)/,, ey (5) 



1,3 = 1 
i?f3 



where e.y is the N x N matrix {eij)ki = 5ik Sji and A G C is the spectral parameter, we find that the 
equations of motion can be written in the Lax form 

(6) 



(7) 
(8) 
(9) 
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when a; — > 0. 





Eq.(H) is the famous functional equation of Calogero. Its general solution is |Q, g 

^(x,A) = -^|±AJ^, Vi.)^pi.) (10) 

where a and p are Weierstrass elliptic functions, the relevant properties of which are recalled in the appendix. 
The elliptic 0{N) Euler-Calogero-Moser model is precisely defined by eq.(^ with V{x) — p{x) together with 
the Poisson brackets (WT 



2 The r-matrix 

From eq.dHI) it follows that trL"{X) is a set of conserved quantities. In particular 



N 



trL{\)=^p,, trL^{\)^2H + p{\). 



The involution property of these quantities trL'^{\) will follow from the existence of an r-matrix which we 
now calculate [g, 0]. Introducing the notations Li{\) = L{X) ® 1 and L2{\) = 1 L(\) we show that the 
Poisson brackets of the Lax matrix elements can be recast as 

{Li(A),L2(/i)} = [ri2(A,Ai),Li(A)] - [r2i(/i, A), LalA^)]. (11) 

Following gl we assume that r is of the form 

N N N 



ri2 (A, ^) = a(A, /^) ^ e^i ® e^i + ^ bij (A, n)eij ® Cji + ^ c^j (A, ^)e 



ij ^/y C-2J . 



!^3 j^i 



Requiring that 7'i2(A, /i) be independent of the pi variables we obtain 

hj{\ii) = -6j.(m,A) (12) 

Cij{X,fi) = Ci-i{fi,X). (13) 

Moreover assuming that ri2(A,/i) is independent of the fij variables yields the following system: 

a(A,^) /(%, A) - bij(X,n) l{qij,p) +Cy(A,^) l{qji,n) = ~l'{qij,X) (14) 

bij{X,^j) l{qjkA) - btk{X,fJ.) l{qjk,^J.) = -2 KlikA) Klii^f^) (15) 

Cij(X,n) l{qjk,X) +Cik{X,n) l{qk],tj) = ^ ^(lik,X) Klij^^) (16) 

Cy (A, ^) l{qki,X) + Ckj{X, /i) Z(g,fc, m) = o ^('^fei' ^) K^tj^f^)- (17) 



A solution to these equations is 



a(A,M) = -i [C(A + Ai) + C(A-M)] 
bij{X,n) = -l{q,j,X~fi) 

C^j{X,^J.) = - l{qij,X + fj.). 



Indeed substituting the preceding expressions in cq.( |l5| , |16|l7 ) leads to the same relation: 



l{qij,X~^i) ?((7j7c, A) +/(qfcj,M- A) l{qjk, tJ.) + l{qtk, >^) l{qji,fJ.)^0 



which upon setting x — ^ (A + qij), y = \ {2fi — A + Qji), z = \ (A + q.j{) and t = hJy^'^ ^ 1>^i + 1>^]) ^^ 
a direct consequence of relation (153). The expression for a(A, /i) is then given by eg. (J14|) , and is simplified 
using eq.M) and (p5[). Finally the r-matrix reads 



N , N 



i,i=l 1,3 = 1 

1 ^ 

- ^ [C(A + n) + C(A - Ai)] 51 ^^^ ^ e,,. (18) 



3 The sl{N) model 

The above 0{N) model can be obtained from the more general sl{N) model by a mean procedure |pl pi Eoj. 
The sl{N) elliptic Euler-Calogero Moser model is defined by the Hamiltonian 

N N 

ff=2Ep'-2E/^^/^'Pfe,) (19) 

i=l i.i=i 

and the Poisson brackets 

fe,g,} = 5,, (20) 

{fij,fkl} = Sjkfil~6iifkj- (21) 

For this model we define a Lax matrix as 

N N 

L{X) =J2{p^~ C{X)h) eu + Y. ^(9«J' ^) /y ey- (22) 

1=1 i,J = l 

The Hamiltonian is given by i7 = ^ / -^^trL^lX). A direct calculation gives 

{ii(A),L2(M)} - [n2{X,^^),L,{X)]-[r2l{pA),L2{^l)] 

N 

- X! ^'(*i' A - A^) (/« - /ji) ey ® e^j (23) 

with the beautifully simple r-matrix 

ri2(A, n) = -C(A - A*) X! "^^ ® Cji + ^ ^(ftj, A - ^) e^j «) ejj. (24) 



i.j=i 



JV AT 



j=l i,j=i 



At this point let us make two remarks: 

• Because of the third term in the right member of eq.(|23|) the integrals of motion trL^{\) are not in 
involution. However we can restrict ourselves to the manifolds {fa = constant)i=i...jv since trL'^{X) 
Poisson-commute with fa. On these manifolds trL'^(X) are in involution. 

• The r-matrix for the 0{N) model eq.(|24|) is immediately seen to be of the form 

where a is the involutive automorphism 

a : A"ey i — > -(-A)"ey. 
This is typical of a mean construction. 



In the following wc will restrict the fij to a symplectic leaf of the Poisson manifold (pi|). Introducing vectors 

{£,z)i=l...N with £,i = {£,°^)a=l-.r 

{r]i)i=i...N with T]i = {■r]1)a=i...r 
with the Poisson brackets 

ur4} = 0' K,'/'} = o= {er,ry^} = -<5y<s.b, (25) 

we parametrize the fij as follows: 

r 

/y—(e»i^,)-E^"'?"- (26) 

a=l 

The phase space now becomes a true symplectic manifold. 



4 The r-matrix of the eUiptic Calogero model 

We show here that the r-matrix for the elliptic Calogero model [g|, Q can be obtained from cq.(p3|) by a 
Hamiltonian reduction procedure [g[ ^ |l^ . 

We choose r = 1 in eq.(g^). On the manifold fij — ^i rjj acts an Abelian Lie group 



Remark that the group acts on L(A) as conjugation by the matrix diag(Ai 



(27) 



■\iji=i---N and therefore all the Hamil- 
tonians trL"'{X) are invariant. Thus one can apply the method of Hamiltonian reduction. The infinitesimal 
generator of this action is 

N 

He — 2_^ ^i fii^ Ai = 1 + ei . 

i=\ 

We fix the momentum by choosing 

fa = OL. 
To compute the reduced Poisson brackets of the Lax matrix, we remark that the matrix 

LC'-'(A) = g-^L(\)g with g = diag(f,).=i-/ 



■N 



N 



N 



^[Pi - aC(A)] eu + a ^ l{qij,\) e. 



(28) 



i=l 






is invariant under the isotropy group Ga of a (which is the whole group itself since it is Abelian) and we 
can compute the Poisson brackets of its matrix elements directly. We find 

{Lf"'(A),L^'^'(A*)} = [r?,^'{X,f,),L?^\X)] - [rg^'{fl,X),L^'^\^,)] (29) 

with 



r^2'^'(A,^)=gr\92- 



1, 



rl2{\^J■) - {guL2i^i)}gl + -[^12,^2(^4)] 



9i92 



where U12 = {5i,<?2}ffi .92 i^ here equal to zero. Redefining 



rf2"' ( A, A*) — ^ rp2"' ( A, Ai) + —J2 e» ® ^u , ^2 (m) 
does not change eq. ([2^) and yields exactly the r-matrix found in |ll|, |lj] 

N N 

r?2'''^{^^ A*) = X! ^(*y ' ^ ^ '-*) ^'3 ® ^i' "'' 9 X! ^(1^l'^^) (' 






N 



■II -r ojjy 






[C(a - /i) -I- C(m)] X! ^" ® ^"• 



(30) 



1=1 



5 Yangian symmetry in the trigonometric case 

The parametrization ( pq ) of fij introduces a sl{r) symnietry into the theory. The transformation 

r 

r 
h=l 

leaves the /ij invariant and therefore also the Haniiltonians. This symmetry is generated by a set of conserved 
currents 

N 

^o' = E^^^ (31) 

It is remarkable that this current was shown, in the rational case W , to be the first of a hierarchy building 
a current algebra commuting with the Hamiltonian — and more generally with a subset of the commuting 
Hamiltonians. 

We now extend this result to the trigonometric case, and we will show that the hierarchy of currents 
form a Yangian symmetry in this case. Taking the trigonometric limit (toi = cx) and (^2 = if) in the above 
formulas, we see that the Lax matrix can be taken of the form 

L(A) = Lo- coth(A)i^ (32) 

with 

N N N 

Lo = ^PiCii - E coth(gy) fij Cij, ^ = E -^^^ ^^^- ^^^^ 

i—l i.j = i i,7 — 1 

By a straightforward calculation, or taking the limit of the elliptic case, we find 

{Li(A),L2(m)} - [r?2,Li(A)]-[r2^,L2(M)] 

-i (1 - coth(A) coth(A.)) ([C, Fi] - [C, F2]) 

^ 1 

- E ^■f" ^ fn)—rj: — rey (g) e^-, (34) 

fj^i smh (q.j) 

where 

N 

r°2 = - E coth((jy) eij ® Cji (35) 






and C is the Casimir element of sl{N) 

N 
C — y. ^ij ® ^ji- (36) 



In spite of the unusual second term in cq.(34) the quantities trL'^{X) are still in involution on the manifolds 
Sq : {fa = a)i=i...N- Indeed, 

{irL"(A),irL™(^)} = „ rn E A__^ [l«-1(a)],^. [L™-i(^)]^, 

f-^i smh (g„-) 



n m 



(1 -coth(A)coth(At)) tri2 (Lr'(A)i™~'(Ai)[C,Fi -F2]) 



and since tr2 ((1 (i$ A) C) — A, we obtain 



i^3 



n m , 



-(1 -coth(A)coth(/^)) ir{L"-i(A)[L"-i(^),i^] -i™-i(/i)[i"-i(A),i^]}. 

If we notice that 

L{\)-L{^i) 



coth(A) — coth(/i) 

we immediately get the involution property. 

We consider now the subset tr{L"-) = tr{Lo + F)" of commuting Hamiltonians; notice that H belongs to 
this subset, since H — \trL^ — a trL + ^Na^ . 

We introduce the following quantities: 

J^^ ^tr{U'F''^), a,6=l,---,r- n = 0, 1, ■ • • , cx) (37) 

where F""^ is the N x N matrix of elements 

(J^"%-/5— eN-. (38) 

We define the generating functional of the currents J^^. It is the r x r matrix T{z) of elements 

T^'iz) = -Is., - Y: ^Jf = -ISa, + tr {j^/''') ■ (39) 

neN ^ ' ^ 

Proposition. On the manifolds Sq we have the foUowing two properties: 

1. The currents J^'' Poisson commute with ah the quantities of the form tr{L'"'). 

2. The generating functional T{z) satisfies the dehning relation of a (classical) Yangian algebra: 

{T{y) ? T{z)} = [R(y, z),T{y) T{z)] (40) 

with 

n '' 

R{y,z) = -2 , n= V eafc^Cba- (41) 

y — z ^-^ 

'^ a, 6=1 

Proof. To prove this proposition we need the Poisson brackets 

^ 1 

{Li,L2} = V\2M\-VlxM\^Y{!^^~!n)—-J--&^®e3i (42) 

tlii smh (qy) 

{ii,Ff} = [-r^i+C,Ff] (43) 

{Fl\F^^} = {5.aFt-5,,Ff)C. (44) 

Remark that the currents Jf^ and the Hamiltonians tr{L") are invariant under the symmetry 

Therefore we can compute their Poisson brackets on the reduced phase space straightforwardly; restricting 
ourselves to the manifolds fa = a, the last term in eq.(E2h vanishes, and we will systematically drop its 
contribution in intermediate calculations. 

We emphasize that in eq. ( p^ , p3| ) the same r-matrix appears. Moreover it is the term [C,i^^''] in eq.(p^) 
which is responsible for the quadratic form of eq.([4G), as we shall see in what follows. 



Introducing the generating functional H{z) = tr{j^) of the Hamiltonians tr[L"') we compute 



^L^—y^^'^T^z^ - 



Taking the trace we obtain 



12 



L2- z L2- z Li -y 



pab 



1 



1 



'21 



pab 



1 



Li-y Li-y L2 - z 



1 1 

Li-y L2- z 



[C,F^ 



a6l 



Lo- Z 



{T'^\y),H{z)}^tr[F-^ 



1 



1 



0. 



L-y'{L~z)\ 

This proves the first part of the proposition. To prove the second part we evaluate 

1 



^ZT^^-Z^^"^ - - 



12 



L2-Z '^ L2-Z 



TT'cd Tpab 

'^2 1 7 — Z^i 



1 



1 



'21 



Tpab 



Li-y 
1 



Li-y Li-y L2 - z 

1 1 



Li - y L2 - z 

1 1 

Li-y L2- z 



{6adFf ~ dcbFf) C 
1 



[C, Ff 



L2- z 



-F 



cd 



Li-y 



Hence taking the trace we get 



{T''''{y),T"\z)} = tr( ^ ^ /j; i,cb 



L — y L — z 



idadF'" - <JchF"'*) 



-tr 



1 



.L-y 

Using the cyclicity of the trace and 



T^cd 



1 



nafo 



L — z L ~ z 



tr 



L- z 



-F' 



ah 



L-y L-y 



, F 



cd 



1 1 



L — y L — z y — z \L — y L — z 



this becomes 



{T'^\y), r\z)] = -^ (SadiT^'iy) - T^\z)) - 6,,{T'^\y) - T'^^z))) 



y- z 
2 



-tr 



Tpcd J^o-^ 



j^ab pct^ 



y — z \L — y L — z L — y L — z 



Remarking that 



tr 



T^cd JT'^^ 



L — y L — z 



N 



L- z 



kl 



i]ki=i ^ ^ ^ «J ^ 



e? v'i 



kl=l "- ^ kl 



b „c 



Ci vl 



T'-'iy) + \5a^ {t^\z) + \5,b 



we prove the result (|40|).D 



The rational limit is obtained by applying the canonical transformation 



Pi 
Qi 



and sending e to zero. In this limit 



io 



'12 



1 

- Pi 

e 

e Qi 



-^rational 

e 

12 rational" 



The Casimir term drops therefore from eq.(43), leaving us with a Hnear Poisson algebra 



{Tiy) f T{z)} = --[Riy, z),T{y) ® 1 + 1 ® T{z)] 



(45) 



which is the result found by Gibbons and Hermsen. 



6 Conclusion 

The Euler-Calogero-Moser model is becoming more and more interesting. On the one hand the computation 
of the classical r-matrix is made considerably easier by the existence of the extra variables fij , the more com- 
plicated r-matrix of the Calogero-Moser model following naturally from a Hamiltonian reduction procedure. 
On the other hand, this model exhibits an exact infinite symmetry which is just a current algebra symmetry 
in the rational case and becomes an exact Yangian symmetry in the trigonometric case. This structure is 
very much reminiscent of the one discovered in [Q . Actually the two currents Jo and Ji (which generate 
the full algebra) are identical in the two cases. Indeed, in our case we have 



N 



N 



jf=T.p^f^~i: 



iir'i smh{q,^ 



£ rah 

-JijJji ■ 



Setting (Xj)'"'' = f^'', e^j = 2^5^-^- and using eq.(|2g) we can rewrite 



N 



N 



jf = Y,p,xf -Y,Q,,(xa,r 



i=l 



i,3 = l 
i#3 



This is exactly the current found in |l3|. In fact the model considered in |lj, |l5|, |l3[ is a quantum version 
of our model for a particular choice of orbit. 

At this point two interesting problems arise. One is the understanding of the role of the r-matrix in the 
quantization of these models. The other is the hypothetical extension of these results to the elliptic case, 
which still remains quite mysterious. 

Acknowledgements We thank D. Bernard for discussions. 



Appendix 



The Weierstrass a function of periods 2ti>i , 2uj2 is the entire function defined by 

2" 



a{z)=z n (l 



exp 



1 



2\uj. 



(46) 



with ujmn = 2muji + 2jiuj2- The functions C and p are 

Ciz) = ^ , Piz) = -C'(^), (47) 

cr(z) 

these functions having the symmetries 

a{-z) = -a{z), C(-^) = -CW, p(-z) = p(z). (48) 

Their behaviour at the neighbourhood of zero is 

a(z) = z + 0(z5) , Ciz) ^ Z-' + Oiz^) , p{z) = z-^ + Oiz^). (49) 

Setting 

%A) = -4i±^ (50) 

a{q) cr(A) 

it is easy to check that 

l{~q,X)^~l{q,~X), l'{q,X) = l{q,X)[C{X + q)-aq)]- (51) 

We need several non trivial relations: 

cr(A - /i) ajX + ^) 
- a^(A) aH^^) = ^(^^ - ^(^)' ('') 

cr(a: — y)(T(a:: + j/)(T(z — t)(T(z + t) + o'(y — z)(T(j/ + z)cr(a; — i)cr(a; + i) +cr(z — 2;)fT(z + a;)(T(y — t)cr(j/ + t) = 0, (53) 
a(2z) a(x + y) a(x — y) ^, ^ ^, n ^, n ,/ x /,,n 

. _^ ; V / 1 X — 7 - a^ +z)- c{x -z) + c(y -z)- ay + ^ , 54 

(T(x + z) a(x — z) (j[y + z) a(y — z) 
this last equation becoming, in terms of the l{q, A) function, 

Choosing the periods uji ~ oo and 0^2 = i? we obtain the hyperbohc case 

6 J ' --v-^-^^^-v^y 3 , ^(^) = — p(^) ■ 3 

and 



a{z) = sinh(z) cxp ( -— ) , ({z) = coth(z) - - , p{z) = ^^^^ ^ ^ + - (56) 



If \\ sinh(A + q) / Xq\ 

smh(A) smh((7) \ 3 ' 



All these formulas were collected in [ pT[ . 
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